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Abstract 



Systems with long-range interactions when quenched into a metastable 

^ ■ state near the pseudo-spinodal exhibit nucleation processes that are quite 

different from the classical nucleation seen near the coexistence curve. In 

^_^ systems with long-range elastic forces the description of the nucleation process 

j-H I can be quite subtle due to the presence of bulk/interface elastic compatibility 

O I constraints. We analyze the nucleation process in a simple 2d model with 

elastic forces and show that the nucleation process generates critical droplets 

with a different structure than the stable phase. This has implications for 

^ . nucleation in many crystal-crystal transitions and the structure of the final 

^ ', state. 
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^ ' Nucleation in systems with long-range forces can be very different ll|-§| than the process 

predicted by the classical theory P,|^. The reason for the difference is the presence of a 
pseudo-spinodal(defined below) p,^] that affects the structure of the critical droplet p|-^ and 
alters the dependence of the nucleation rate on the thermodynamic parameters ||l|;10]. An 
interesting and important class of materials that exhibits pseudo-spinodal behavior is those 

^ . that interact through elastic forces; a subclass of which undergoes martensitic structural 

Q I transitions. []TT| An example of such a transition is when alloys such as FePd and NiTi 

O 



transform on cooling from an "austenite" phase at high temperatures to an equal width 
mesoscale twin phase below the martensite transition temperature To |T2|. This transition 



^ ■ is first order and takes place via nucleation. However, the nucleation process in these systems 

j^ . is not well understood [|T^-|T^. 

The purpose of this Letter is to present an analysis of nucleation near the pseudo-spinodal 
of a model with elastic forces. A complete specification of the critical droplet profile requires, 
in addition to the usual bulk terms, [0,|1O| a consideration of a term generated by the elastic 
compatibility constraints. This is, to our knowledge, the first treatment that describes 
nucleation near the pseudo-spinodal in this class of materials and is the first indication that 
the compatibility constraint plays an essential role in this nucleation process. In addition, 
this is the first indication that homogeneous nucleation in a crystal- crystal phase transition 
takes place with a critical droplet that does not have the symmetry of the stable phase. 



In this Letter we analyze a model that exhibits a two-dimensional square to rectangle 
transition [16| and captures the essential physics of crystal-crystal transitions in systems 
with elastic forces. The order parameter(OP) is a rectangular or deviatoric strain, which is 
a symmetry adapted combination of the 2d strain tensor £:^j,(/i, v = x,y). The non-OP, or 
secondary, strain components are related to the OP through a compatibility equation. We 
can write a Ginzburg-Landau (GL) free energy in the form |jl6| F(e) = -Fo(e) + -Fc,rad(Ve) + 
-Fes (61,62), where 

F„(e) = J df\{T - l)e\r) + e\r){e\r) - 1)% (1) 



grad 



(V6) = /rff[^(Ve(f))2 + -{V\{r-)n (2) 



FUei,e2) = ^Jdf[A,el{r) + A2el{f)]. (3) 

Here e(f) = -^[£{xx)ir) - £(yy){r)] is the OP strain, ei(r) = ■^[£{xx){r) +i£{yy){r)], e2(r) = 
£(xy)i'i^ are the compression and shear strain, respectively and r = jf—jfr, where T^. is the 
temperature at which the OP would completely soften. The subscripts (xx), (xy) and (yy) 
denote partial derivatives, Ai and A2 are elastic constants for the compression and shear, 
and a and b are strain gradient constants independant of T. 

The St. Venant compatibility equation for the symmetric strain tensor is V x (V x 
£{r))'^ = 0. Using the Lagrangian multiplier formalism W^ in d = 2 we find for Fourier 
expandable strains e{k) that ei{k) and 62 (A;) are proportional to e{k), the Fourier transform 
of e(r) with k dependent coefficients. This result allows us to replace Fcs{ei,e2) in Eq.(^ 
with an OP potential Fc,(e) = F^f''{e) + F^7^^"^'=(e) where MM 

Flf\t)= [ dkU'"''\k)\t{k)\', (4) 
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U'-'\k) = ^"^ ^\J, , (5) 



where the surface term is generated by the interface between the high symmetry "austen- 
ite" phase and the low symmetry "martensite" phase. In Eq.(|^) the integral is over the 
momentum modes of the interface, J is independent of k and depends only on the ratio 

p = A2/2A1, and 

^ ~^ Pi2\kxky\\ 



In the above equations ^^.(y) is the x{y) component of k with respect to the underlying square 
lattice |18[ and k = \k\. 



In the thermodynamic limit we can ignore F^^^f"-'^'^ and we first assume a spatially ho- 
mogeneous e(r'). The free energy has the following form; For r > 4/3 there is one minimum 
at e = 0. For 1 < r < 4/3 there are three minima but e = is the global minimum. 
For < r < 1 there are again three minima. The one at e = (the "austenite" phase) 
is metastable while the two minima, symmetric about e = (the "martensite" phase) are 
stable. For r < there is no longer a minimum at e = 0. Note that this analysis does not 
consider instabilities to perturbations with non-zero wave vectors. We will return to this 
point below. We will first investigate the nucleation process from the metastable minimum 
at e = near the spinodal at r = 0. From the spatially homogeneous free energy it is 
simple to calculate the order parameter exponent e ~ r^ with f3 = 1/2. Reinstating the 
Laplacian term in free energy it can be seen that the correlation length diverges as ^ ~ t'" 
with u = 1/2. Adding a spatially homogeneous external field to the free energy leads to a 
susceptibility x ~ 't~'^ with 7 = 1 [p^ . 

Turning to nucleation near r = 0, we first note that nucleation does not occur in systems 



with infinite range interactions. That includes mean-field systems |T0|. However, systems 
with elastic forces do not have infinite range interactions due to the screening from defects. 
pO| Hence we will be dealing with systems that have the same bulk and surface interactions 
as in eqs.(||-|]) except they will have an exponential cutoff of the form exp(— -^), where 
i? >> 1 is the interaction range and r = |r|. 

Since the interaction range is large but finite the system is no longer mean-field but near- 
mean-field [0,^ • In order for the mean-field approach, including the idea of a spinodal, to 
be a reasonable approximation when R ^ 00 the system must satisfy the Ginzburg criterion 
2l|,|22| namely 



M-2 = —-r- « 1- (8) 

Note that the correlation length,^, as are all lengths, is in units of R. The Ginzburg criterion 
can be rewritten as R'^r'^-'^/'^ = A >> 1 where A is a fixed large number. When the Ginzburg 
criterion is satisfied, many aspects of the mean-field spinodal are still present. However, the 
singularity has been smeared out. |||,P] The larger A the better the spinodal is approximated 
by the pseudo-spinodal. Since A » 1 for these systems the pseudo-spinodal is very close 
to a true spinodal. 

To calculate the nucleation or critical droplet structure we will use saddle point techniques 
il|J^J^,p!0[] . Near the pseudo-spinodal there will be an incompletely softened mode that can 
be identified by examining the k coefficient of the Gaussian term in the action. We take the 
action to be the free energy in Eqs.(||-|^) ^. Initially we will ignore the surface term given 
in Eq.(|). The structure factor S{k) is then 

S(k) ~ [r + ^e + h^ + U'^'\k)]-\ (9) 
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— * 
where the +(-) is for a positive(negative). Consider first a, 5 > 0. Since all terms in S{k) 

are positive semi-definite the only divergence is when fc —* and k^ = ky. The surface 



term in Eq.(|^) would appear to strongly suppress the fluctuations that cause the divergence. 



However, this term was derived for a sharp interface [[16ip!8|. We need to extend this result 
to an interface with width ^. The reason for this particular scale will become clear. We 
can consider the smooth interface to be a sequence of sharp interfaces or steps each one 
contributing a term to the free energy of the form given in Eq.(^. Since k^ = ky, I{k) = 0. 
We will assume a form for the interface of emterface oc exp(— r/^) and C, » R. The difference 
in order parameter amplitude between steps is 

(i exp(— r/^) 
Aeinterface oc — exp(-r/4)dr = R, (10) 

where we have chosen our differential step dr = R. This is justified since all lengths are 
in units of R and R is the coarse graining scale in the GL theory @,|10|. Taking k ~ ^~^ 
and dk ~ ^~^ in Eq.(^ the free energy cost for the surface is F^'^^f"-^^ ~ ^ g , where ^e{r) 

scales as e{k) and the number of steps in the surface is C,/R- As we will see the dominant 
contribution to the free energy F scales as J dfTe'^{r) = ^'^re'^, where by e we mean to 
include only the dependence of e{x) on r. Comparing these two scaling forms we have 
psurface _ F^j^ << F, whcre C is a constant of order 1 and A » 1 from the Ginzburg 
criterion. The Ginzburg criterion also implies that for a fixed but large R, r cannot reach 
zero 1 10] and the system remain near- mean- field. The ratio p^^'^f'^'^^^ j p = Q only as R and 
hence A ^ oo. For the long-range potential we are using the surface term is, in general, 
small enough to neglect. Note that if A is not infinite then the surface term is added to S{k) 
eliminating the divergence at A; = 0. Since A » 1 the structure factor can be extremely 
large and the true spinodal is well approximated. In calculating the surface contribution we 
have assumed that there is a domain with non-zero e(x) with a linear size of the correlation 
length C, imbedded in the metastable e(x) = phase. We now proceed to demonstrate the 
existence of this domain. First we note that a and b in Eq.(^ must have units of length to 
the second and fourth powers respectively and hence are proportional to R^ and R'^ as all 
lengths must be proportional to R ||10[- The Euler-Lagrange equation for the critical droplet 
is obtained by setting the functional derivative of F(e) in Eqs.(|}]^) equal to zero to obtain 

- -Vh{r) + ^V^e(f) + 2re(f) - Ae^r) + 6e^(f) + / rfr=^t/^"''=(f, r^)e(r=') = (11) 

where [/(r) is f/^""^(r) multiplied by the exponential cutoff exp(— -^). We now assume a 
solution of the form 



e(0 = Ecn(r)exp(z£,„ ■ f)V^(y) = G(r,r)7/;(^), (12) 

L/ Li 



where L >> R, koQ is the value of k^^n at which the mean- field structure factor (Eq. 
diverges and co(r) >> Cn{T) for n 7^ 0. For r ~ the c„ for n 7^ can be neglected p2| . 
We are near a pronounced pseudo-spinodal so that we expect the critical droplet to have 



an interior structure similar to spinodal critical fluctuations [|TOl(See Eq.(||).). Since a > 
implies kofl = 0, with the assumed form for e(r) in Eq.(0) the Euler-Lagrange equation 
becomes 



- f VV(^) + ^VV(^) + 2rV^(^) - Ai;\^) + 6^'{^) = (13) 

Since kx = ky the term involving f/''""^ gives no contribution and the Cnir) are chosen so 
that G3^r,r) =G{T,f). 

Since r ~ the solution of Eq.(O) has the scaled form 



V;(f) = Dr'/^A (14) 

where the V'^^(^) and ip^ij^) terms have been neglected since they are higher order in r, 
B and D are constants that can be determined from Eq.(|T3|) and L = ^. This form of 
the solution is what we assumed when we calculated the contribution of the surface term 
in Eq.(§) for a smooth interface. Hence the omission of the surface term is justified self 
consistently as is the scaling of the bulk free energy used to compare with the surface 
contribution. 

The nucleation barrier, AF, is calculated by inserting the critical droplet solution, 
Eq.(|l2]), into the free energy 0, Eqs.(|ll-|D. It is straightforward to see that AF oc Rd.^2-d/2 ^ 
A. Therefore ii A = oo, the system is mean-field, rather than near- mean-field, and there is 
no nucleation. 

Note that the saddle point object which is the nucleation droplet shows no evidence of 
the twin stripes seen in the simulation of the stable phase of this model [0. The critical 
droplets near the pseudo-spinodal are unstable p|,p!0| and differ from the metastable phase 



by an order of magnitude given by r^/^ ~ 0. Their initial growth phase is a "filling in" or 
an increase in the order parameter difference. p|,p!0| The filled in droplet will have a sharp 



interface and hence must have twinning |TB|]. Therefore, the symmetry breaking which results 



in the twin stripes must appear in the growth phase. As we will see, if a, the coefficient of 
the (Ve(r))^ term, is negative the case is somewhat different. We treat this next. 

For a < and 6 > we take the minus sign in the the structure factor in Eq.(|^). Since 
jjbuik^f^-^ is independent of k, the value of k where the structure factor diverges is 

fc = ±i[M±(H!_(L±i^)i/2]i/2 (^5) 

2U M6 2 ^ ^ ' ^ ^ 
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where < Tq < 4p is a fixed value of U^^^^{k). For r > ^-^ — r^ there is no divergence for 
real k and hence no instability. Since the largest value of r for which there is an instability 
is the spinodal then for a < 0, the spinodal is at r^ = -^ > 0. The structure factor will now 
diverge at a non-zero value oi k = kg where ko is given by Eq.(^) with r^ = 0. Note that 
the additional instability generated by a < is at a value of r greater than r = expected 
from a simple thermodynamics calculation. It is straightforward to calculate the exponents 
of the correlation length and the order parameter which have the same values as those at 
the r = spinodal for a > 0. 

Turning to the nucleation problem for a < and initially ignoring the surface term, the 
Euler-Lagrange equation has the form 

/ df^S-\f^)e{f- r^) - Ae\r) + 6e\r) = 0, (16) 
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where S~ (r) is given by the Fourier transform of the inverse of S{k{) in Eq.(|^) with U (k) 
replaced by U^^^^{k). We again assume a solution of the form given in Eq.(^) where kofl is 
the vector at which the structure factor diverges with a < 0. Since Tq = implies k^ = ky and 
l^o.ol is given by Eg. ([T5|) , kofi is specified. We now expand i^i^-^) in a gradient expansion 
about ^ to obtain 

^12'7'^i^) + (r - r,)7o^(^) - 4^'(^) = 0, (17) 

where, anticipating the scaling, higher order derivatives and higher powers of '?/'(^) have 
been neglected. The constants 70 and 72 are J dr[exp{ikofl ■ r) + exp(— -ifco^o ■ r)]S^^{r) and 
/ drr'^[exY){ikofl ■ r) + exp(— ifco ■ ^]S~^{r)., respectively. As above, the c„ have been chosen 
so that G^{t^ r) = G{t, r). The solution of Eq.(17) is of the form 



^{-) = A{r - r^y/'^i—) (18) 

justifying the omission of higher order terms. 
The strain field of the critical droplet 

e{r} ~ A{t - T,)^/^[exp{tK,o ■ ^ + exp{-tK,o ■ r)]V^(-^), (19) 

where we have neglected terms with n 7^ 0, is not that of the stable phase but does exhibit 
a spatial modulation of regions where e{r) ^ 0. Consequently, the stable phase structure, 
as in the case a > where there is no spatial variation in the strain, must evolve during the 
growth phase. Note that the solution, Eq . ([T8|) , justifies the omission of the surface term via 
an argument virtually identical to the one given above. 

We have calculated the first critical droplet structures for nucleation from an "austenite" 
like phase to a twinned "martensite" like phase near the pseudo-spinodal in a system with 
elastic forces. The droplets do not have the stable phase structure as expected from classical 
nucleation [^0 and in the a > case exhibit no spatial modulation. Droplets that do not 
have the stable phase structure have been predicted in the nucleation of the crystal from 
the melt [§,^ but this is the first indication of such a droplet structure in a crystal- crystal 
transition. It is also the first result that demonstrates the importance of the compatibility 
constraints to the phase transition kinetics. 

It is important to note in systems with i? >> 1 that classical nucleation is strongly 
suppressed. In the classical case the nucleation rate is proportional to exp(— ^Jli ) B|31 
where a is the surface tension between the droplet and the surrounding metastable state 
and A/ is the free energy density difference between the stable and metastable states. For 
classical nucleation, near the coexistence curve, a ~ 1 so that for i? >> 1 nucleation is 
severely suppressed. In order to have nucleation in a reasonable time frame the quench 
must bring the system close to the pseudo-spinodal where a << 1. Therefore, nucleation 
near the pseudo-spinodal will dominate the phase transition process in realistic experiments. 
Finally we note that the form of nucleation discussed in this Letter allows the possibility of 
evolution into metastable crystallites with symmetries different than the stable phase. 
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